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Abstract: Studies of models of current flow behaviour in Electrical Impedance Tomography 
(EIT) have shown that the current density distribution varies extremely rapidly near the edge of 
the electrodes used in the technique. This behaviour imposes severe restrictions on the numerical 
techniques used in image reconstruction algorithms. In this paper we have considered a simple 
two dimensional case and we have shown how the theory of end point / pinch singularities which 
was developed for studying the anomalous thresholds encountered in elementary particle physics 
can be used to give a complete description of the analytic structure of the current density near 
to the edge of the electrodes. As a byproduct of this study it was possible to give a complete 
description of the Riemann sheet manifold of the eigenfunctions of the logarithmic kernel. These 
methods can be readily extended to other weakly singular kernels. 



I. INTRODUCTION 

There are numerous examples of practical situations where electric current is used to probe 
the interior of some object of interest. One emerging technology which specifically uses this 
approach has become known as Electrical Impedance Tomography (EIT). This is a method of 
medical and industrial imaging in which electrical currents are applied to the surface of an object 
and the induced surface voltage is measured. These data are then used to produce an image of 
the conductivity distribution in the interior of the object. An extensive literature exists on EIT.^ 

The particular feature of EIT which is of interest to us here is related to the observation 
that in practice the electric current can be applied only through a finite number of electrodes — 
currently in the range 16-64 for two dimensional applications. The consequences of this fact and 
the appropriate mathematical modelling of the electrodes have been discussed in a number of 
papers. 

In medical applications one of the significant problems for EIT is the existence of a thin layer 
of material of high, but unknown, contact impedance lying between the current drive electrodes 
and the body. Various models have been proposed to describe this phenomenon but one which 
has strong experimental support^ is to suppose that on the (finite size) electrodes the electric 
potential, $, is related to the electric current a{d^/dn) by 

$ + Za— = V 

on 

where a is the conductivity just below the electrode, Z is the contact impedance and V is the 
potential of the electrode (a constant). The induced voltage, $, is found by making measurements 
on high impedance electrodes, also attached to the surface. 

This model has been studied numerically using a boundary Fourier technique^ for the case 
of constant Z and an interesting phenomenon which was observed was the appearance of very 
sharp peaks in the current density distribution at the edge of the electrodes. More recently 
the model has been studied numerically for non-constant Z using the weakly singular integral 
equation described in Eq. (j^) below. Although the singularity is weak, its existence has important 
consequences for the numerical treatment of this equation. Details of this work will be given 
elsewhere^ but the point which we wish to emphasize here is that the distinctive sharp peak 
behaviour occurs for a wide range of conductivity distributions and contact impedance forms. In 
Fig. 1 we show typical results for the current density distribution in the case of eight electrodes 
with an input current on the Ith electrode of cos 6i. 
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FIG. 1. The current density distribution obtained solving Eq. (^) numerically for the case 
of 8 electrodes, each having a contact impedance equal to 0.22. 

An unwelcome consequence of the sharpness of these peaks is that the direct numerical 
modelling of the potential with the finite element method has become an excessively substantial 
task due to the high number of mesh points needed near to the edge of the electrodes in order 
to accommodate the rapid variation of the normal derivative of the potential. The aim of this 
article is to give an explicit analytic description of these singularities in a form which should 
substantially improve the speed of the numerical computation. 

Since the appearance of the peaks is a boundary phenomenon which is very little influenced 
by the actual values of the conductivity a inside the disk — see the discussion in Section III about 
the dominant singular integral equation — we shall focus our attention on the constant cr case. 
In this case the governing partial differential equation is Laplace's equation and our problem 
becomes one in potential theory. Consequently, we shall study the electrode model defined above 
for the standard domain of the unit disk. The importance of the unit disk stems from the fact 
that the potential problem for any simply connected two dimensional domain can be reduced to 
the unit disk by an appropriate conformal mapping. 

In earlier investigations of these peaks it had been shown^ that this problem can be solved 
explicitly in the zero contact impedance case = 0), for the case of two electrodes. In this 
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specific case one finds that near the edge of the electrodes 

a$ 1 

dn ^Jx 

where x is the distance, along the boundary, from the edge of the electrode. Thus the normal 
derivative of $ becomes infinite at the edge of the electrodes. However if the contact impedance 
Z is not zero, although d^jdn still has sharp peaks at the edge of the electrodes (see Fig. 1), it 
remains finite since both $ and V are finite and since d^/dn= (V — ^)/{Za). This shows that 
the nature of the singularities in the Z ^ Q case cannot be obtained from that encountered in 
the soluble Z — model. Thus, in order to have a correct understanding of these singularities 
some deeper investigations are necessary and this represents the main goal of the present paper. 

Overview of the paper 

After a short description of the mathematical model in Section II, in Section III we formulate 
the boundary problem as an integral equation. As shown there the kernel of this integral equation 
has a weak (logarithmic) singularity, and this has direct consequences for the singularities of the 
solution near the electrodes edges. We will describe these singularities in terms of an asymptotic 
series for the potential. 

In pursuing this program we shall have to step off the real axis into the neighbouring com- 
plex plane. This will be necessary since we will write the solution of the integral equation as an 
infinite sum of the free term and the eigenfunctions of the weakly singular kernel. As one knows, 
there are many examples of series uniformly and absolutely convergent on the real axis which 
cannot be differentiated term by term [for example cos(nx)/n^ converges uniformly, while 
the kih. derivatives of its terms contain factors of the form n^~'^ which spoil any convergence] , but 
in deriving asymptotic expressions we will often have to perform this kind of operation. However, 
in contrast to what happens on a real interval, in the complex plane there exists a marvellous 
theorem which states that given a sequence {/j} of functions holomorphic in some domain Vt 
which converge uniformly, — > /, on all compact subsets of then (a) / is a holomorphic 
function in Q and (b), // as well as the higher derivatives //"^ tend uniformly towards /' and 
/*^") on any compact subset of Jl. In some way by walking in the complex plane around the 
singularity one has a better view of what really happens there. 

In studying the analytic properties of the free term (Section IV) and of the eigenfunctions 
(Section V), we shall use techniques similar to those from the 'pinch and end point singularities 
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theory' which was developed some time ago by Eden, Landshoff, Ohve and Polkinghorne^ in 
elementary particle physics. However our problem is more complex than that related to the 
Feynman graphs in two respects. First we will have to consider moving cuts rather than moving 
poles, and second, we will no longer have integrals over some explicitly given functions, but 
integrals over the a priori unknown eigenf unctions whose singularities we are trying to find. 

Handling infinite series can also be dangerous because spurious singularities may creep in 
as happens, for example, with the common geometric series. The proof that this does not happen 
in the neighbourhood of the singular point is given in the first subsection of the Section V. The 
analytic properties of the eigenfunctions and the recursive procedure to compute the coefficients 
in the asymptotic expansions are given in Section V B [see Eq. (|6^)]. From this expansion it 
follows that the spikes of the current density near the edge of the electrodes are of the form 

^ (^rnkx"" ■ \og\x) + 0{x''+^-') + regular part, 

'^'^ m=l k=l 

where x is the distance along the boundary from the edge of the electrode and Cmk some real 
coefficients. Since the derivative of xlog(a;) is 1 + log(x), this means that although these spikes 
are finite they have infinite derivatives. 

II. THE MATHEMATICAL MODEL FOR EIT 

The usual model used to describe the forward problem in EIT is obtained by considering 
the object as consisting of isotropic material with conductivity distribution a contained in an 
open, simply connected region Q surrounded by a reasonably smooth boundary dQ. On the 
surface, dfl, a number, L, of electrodes are attached and electrical current is apphed. 

In this case Maxwell's equations give: 

V ■ (ctV$) =0 in n. (1) 

Further, the total current driven on the Ith electrode J; = /p^ cr9$ /dn is a known quantity and 
there is no current outfiow outside the region covered by the electrodes, F = Fi U r2 U .... U F^. 
If we now introduce our electrode model mentioned earlier to the case when a is constant and Q 
is the unit disk, the physical problem is equivalent to the mathematical problem of solving the 
following boundary value problem: 

V^$(z) =0 in (2) 
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-\Vi-^{z)] on the electrode C 9ri, 



(3) 



for constant induced voltage, V/, and total current driven, J^, on each electrode. Here Zi represents 
the contact (the 'skin') impedance and $(2; = e*^ G Ti) is the potential just underneath 'the skin'. 

III. THE INTEGRAL EQUATION 

If the values of the normal derivative 

dn 

were known everywhere on the boundary dVl of the unit disk, we would be considering a classical 
Neumann problem, which is readily solved by means of the formula 



on 



d9' + const. 



(4) 
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where the Neumann kernel N'{z, z'] 0) is 

7V(z,/;0) : 

with z' = e*^' on the unit circle. 

If we integrate the kernel with the values of the normal derivative d^/dn on the unit circle, 
we obtain a function $ which is harmonic throughout the unit disk, which vanishes at 2; = and 
which has the prescribed normal derivative values on the boundary. Similar kernels ^/{z, z'; zq) 
producing functions vanishing at z = zq rather that at z = can be written easily,^ but the 
knowledge of the normal derivative determines $ only up to a constant. 

For what follows it is interesting to continue the Neumann integral (|^) up to the boundary. 

Since 

e^^-e^^'p = 2[l-cos(^'-^)], 
if z and z' are of the form z = e*^ and z' = e'^' the Neumann kernel on the unit circle reads: 



Ar(e^^e^^';0) = --1 log{2[l - cos(0' - ^^)]}, 

ZTT 



log 

TT 



2 sin 



'9' -9' 



(5) 
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In our case the normal derivative is known explicitly only on that part of the boundary which 
lies between the electrodes (i.e., on 9n\r) where d^/dn is identical to zero. However, since the 
integral representation (^) can be continued up to the boundary, conditions @ yield a linear 
integral equation for the boundary values p{6) = $(e*^) of the potential: 



1 L „ jn/ 



Although the kernel 



2 sin 
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p{9'). (6) 



log 



TC 



2 sin 



'9' -9' 



becomes infinite each time 9' equals 9, this logarithmic singularity is weak enough to be 
integrable. The kernel is therefore of Hilbert-Schmidt type and so one can benefit from all 
the advantages of Fredholm integral equations of the second kind, namely the existence and 
uniqueness of an solution p{9). 

Equation may be rewritten in a form which exhibits the logarithmic singularities of the 
kernel. Taking e*^ G Ti^ we may write 

. f9'-9^ 



piO) 



1 



E / log 



2 sin 



'9' -9^ 



[Vi - PjO')] 
Z,(9') 



d9' 



1 

TT Jr,, 



log 



sm 



9' -9 



d9' 



(7) 



In order to study the local behaviour of the solution p{9) near the edges of the l^ih electrode, 
we absorb the first two terms which are continuous into the 'free term' of the so called dominant 
singular integral equation,^ which is of the form: 



f{x)=g{x) + A / Kix,t)fit)dt, 
Jo 

K{x,t) = \og\t 



< a; < 1, 



x\ 



(9) 



Here the points a; = and x = 1 correspond to the edges of the electrode under consideration. 
Equation (||) can be readily deduced from (|^ by a suitable change of variables and functions 
in the Zi{9) constant case, but as we shall show elsewhere^° the discussion for the general case 
(non-constant Zi, non-constant a) is fairly similar. When we follow this procedure we find that, 
as well as the first two terms from Eq. (0) which are regular, the function g{x) contains the term 



TT 



d9' 



Z,. 



-log|^'-^| 



so that after the changes of variables g{x) has the form 

g(x) = / log It — x\w(t)dt + regular part. 
Jo 

For the convenience of some subsequent proofs we shall also be interested in the iterated equations 
obtained by replacing f(t) under the integral by the right hand side of the integral equation (||): 



, 1 

2 



fix) = g2ix) + X' / K2{x,t)f{t)dt 







fix) = g3ix) + X' / Ksix,t)f{t)dt 



(10) 



and so on, where 







g2{x) = g{x) + \[ K{x,t)g{t)dt, 
Jo 



,»i (11) 
gaix) = g2{x) + \^ K2ix,t)g{t)dt, ... 







and 



K2{x,t) = f K{x,T)K{T,t)dT. 

Jo 

K3{x,t) = f K2{x,T)K{T,t)dT. 

Jo 



1 (12) 



If this iteration had been continued indefinitely we would have found the Neumann series for /. 
Since these series usually converge only for very small values of A, we shall not use them but stop 
after a finite number of terms and take advantage of the fact that the eigenvalues of Kj {x, t) are 
the powers {A{} of the eigenvalues {A„} of K{x,t). Indeed this will be quite helpful in some 
subsequent convergence proofs. 

IV. SINGULARITY OF THE FREE TERM 

In this and in the next section we shall try to find the analytic structure of the edge 
singularities of the solution without solving the integral equation, the latter being possible only 
numerically or in some very special cases. ^ To this aim we shall use methods similar to those 
from the theory of the pinch or of the end point singularities/ well known to particle physicists 
working in analytic S-matrix theory. As a preparation to what follows it is probably helpful 
to look to the corresponding chapters from the classical book of Eden, Landshoff, Olive and 
Polkinghorne.^ 
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In this section we shall deal with the free term g{x) of Eq. (j^). As mentioned in the 
Introduction, in order to find the analytic structure of the singularities we have to step into the 
neighbouring complex plane. We shall start our investigations with some negative values zq of 
z for which log |t — -z] = log(t — z) since the integration variable t is between and 1. Having 
to perform analytic continuations we prefer to handle holomorphic expressions (log It — z\ is not 
a holomorphic function of z) and so, instead working with g{x) we shall focus our attention on 
functions of the kind 

F{z) = \og(t - z)w{t)dt (13) 

JO 

where the weight w{t) is a function of t, which is holomorphic (no cuts or other singularities) in 
neighbourhoods of t = and t = 1. Although the holomorphic extension F{z) is different from 
g{x), it is closely related to it since, up to regular terms. Re F{x ± ie) = g{x) for x G [0, 1] and 
e\0. Since the weight w{t) may differ very much from one case to another, integral ( |T3|) cannot 
be performed explicitly. Therefore it will be interesting to have a mathematical procedure which 
should be able to predict the form of the singularity without actually performing the integrals. 
However, in order to have a partial check of the results which will be obtained below, we note 
that in the simplest case w = 1 we obtain Fw=i{z) = z\og{—z) + {1 ~ z) log(l — z) — 1. 

As mentioned at the end of the last section, if A is small enough the Neumann series converge 
and so the solution can be written in terms of the free terms gj of the iterated equations. Therefore 
at the end of this section we shall discuss briefly the singularity of the iterated functions gj{x) 
since they provide a check of the results obtained in Section V in the general case. 

A. Different ways of defining a cut 

We recall that the features which are important when considering the cut structure of 
complex functions are the locations of the branch points and not the way in which the cut is 
taken. Indeed, the cut can be deformed or trailed as will become apparent below, in Figs. 3 to 
5. As a first example, consider the function F^=i{z) given above. To define the cut of the first 
term zlog{—z), we first introduce the function Z{z) = —z and ask then that the cut of log Z{z) 
should run along the positive Z real axis. One achieves this by writing Z = \Z\exp{i(j)) and 
requiring (p to be in the range [0, 27r). With these conditions the cut of the first term, zlog{—z), 
of -^«,=i(-2) runs along the negative z real axis. The function z\og{—z) is real and equal to 
zlogj— z| = — |z|log|z| below the cut. Above the cut it will contain an additional imaginary 
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part equal to 2iT:z. There will be no cut along the positive z real axis, here the value of z log(— 2;) 
is equal to 2; log \ z\ + mz both above and below the real axis. 

It is not compulsory to take the cut of the second term of Fw=i{z) to the left. Indeed, if we 
define a new variable Z[z) — 1 — z. wc can redefine the 'fundamental Ricmann sheet' of log Z by 
requiring that the argument of Z to be between — tt and tt rather than between and 27r, so 
that the cut of logZ — log(l — z) will run along the negative real Z axis. Summarising Fyj=i{z) 
has no cut on the interval [0, 1] but two cuts running from —00 to and from 1 to 00. 
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Fig. 2b. 

FIGS. 2. Example of two different definitions of the cut of log(2; + 1). 

The possible patterns for the cut of the logarithm are not exhausted by the cases discussed 
above. To have a further example let us consider for instance the function Z{z) = 1 + z and 
begin with a cut of log Z running along the negative Z real axis. In this case the cut of log(2; + 1) 
looks as in Fig. 2a where the ±i7r zones mean that the values just above or beneath the cut 
differ by ±i7r from the mean value across the cut. (For illustration purposes we have slightly 
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deplaced the cut: as it stands it corresponds to the function log [2; — (—1 — ie)], £ > 0.) 

However we could alternatively define the cuts to run as in Fig. 2b where again the 
iivr zones mean that the values of the logarithm differ by iivr from the mean value across the 
cut (which is not necessarily real). With these specifications one finds immediately that the value 
of log{z + 1) at the points A, B and C are, up to epsilons, equal respectively to log 1^;^ + 1|, to 
log |za + 1| + 2z7r and again to log |za + 1|- As a further example, the values of log(z + 1) at the 
points D{z = — i/2 + ie) and E{z = — i/2 — ie) are respectively log y^5/4 — (26.5/180)i7r + 2i7c 
and log ^5/4 — (26.5/180)z7r, while the mean value across the cut is there equal to log y^5/4 — 
(26.5/180)i7r + OT. 

B. 'Hooking' the integration contour 

Let us come back to the function F{z). In what follows it is important to consider sepa- 
rately the parameter (the 'control') z-plane, and the complex t-plane where the integration is 
performed. As stated above, our aim is to find the singularities of F{z), in the 'control' 2;-plane, 
without performing the integration explicitly. Since our final goal is the description of the sin- 
gularity of F{z) at 2; = 0, in this section we shall consider only analytic continuations performed 
in some neighbourhood of the origin. A similar discussion can be made for the other end point 
z = l. 
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FIG. 3. Integration path (full line) and cut of the logarithm (dashed). 

Suppose that initially z lies somewhere immediately below the negative real axis in the 
complex z-plane: z = zq = xq — ie, where xq < and e > 0. The function log(t — z) from the 
integral of Eq. (|l3|) has, as a function oft (i.e. in the complex t-plane where the integration is 
performed), a cut running parallel to the real t-axis from t = —00 to t = z (see the dashed line 
in Fig. 3). From the point of view of the integration t-space, z is a parameter. Suppose now that 
z moves towards a point Za = Xa — ie {xa > 0) and then returns to zq {xq < 0) without crossing 
the integration contour. Correspondingly, the head of the cut of the logarithm (as a function 
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of t) will move in the i-plane as shown in Fig. 3, but, since it will never cross the integration 
contour, the value of F{zo) will be identical with the value that the function F{z) had before 
the point z began to move from zq- 
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FIG. 4. The integrals over [0, 1] and over the deformed path are identical if w{t) is holomorphic. 

The situation is however different if the path followed by the point z in the control complex 
z-plane crosses the segment [0, 1] before returning to the initial position z — zq. Here again the 
path followed by the 'head' of the cut begins at Zq and ends at the same point, but this time it 
winds around the integration end point t — 0, crossing the real axis ai t — zi, as shown in Figs. 
4 and 5. 



- inzone 

-Hnzone ^ Zp -unzone y 


Imt 

^^ ^ 


- inzone - inzone 


Ret 



FIG. 5. The moving cut deforms the integration contour and produces an Anomalous Threshold. 

A genuine analytic continuation of the function F{z) should, of course, be at least contin- 
uous, i.e. have no jumps or other discontinuities. Therefore, even before the end point t — z of 
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the singularity of the logarithm reaches the point t = Zi, which lies on the real axis between and 
1, (see Fig. 4), we shall use the freedom we have to deform the integration contours inside the 
analyticity domain of the integrand without altering in any way the value of the integral F{z). 

Since the branch point t = z 'trails' behind it the cut of the logarithm when z moves 
further through points z^ in the upper half 2;-plane (Fig. 4) towards zq, the value we obtain for 
the analytic continuation F^^\zo) of the initial integral will be given by the integral along the 
path C in Fig. 5 (the dashed line in Fig. 5 represents the cut of the logarithm; as it stands, the 
cut ends at the conjugate point zq = xo + is rather than at zq). Since the value of the logarithm 
on the lower lip of the emerging part of the cut contains an additional 2Tii with respect to that 
on the upper lip (recall the discussions about the points A and B from Fig. 2b), the integral on 
the part of the contour around the cut (see Fig. 5) is 

/ {log \t — zq\ + 2i'K)w{t) dt + / log \t — zo\w{t)dt = —2iTT / w{t)dt. 

Jo J Zl) J ZQ 

Here we have supposed that the weight w{t) has no singularities at 2; = 0. Such a point Zq where 
the integration starts and which is below the initial integration threshold is called an 'anomalous 
threshold' in elementary particle physics. 

The new value F'^^\zq) of F{z) obtained by means of this analytic continuation process is 

hence 

F^^\zq) = 2m r w{t)dt + \og{t - Zo)w{t)dt (14) 
JO Jo 

where the logarithm has the same determination as in Eq. (|T3|) (i.e. a cut like in Fig. 3). 

There are many ways of defining the Riemann sheets of F{z) and here we describe only 
two of them: 

(i) If we now define F{x + iy) for any point in the upper half 2;-complex plane to coincide 
with the function F^-^\x + iy) defined above, i.e. F{x + iy) = F^-^\x + iy) (y > 0), we have 
implicitly required that the function F{z) should have no cuts on the segment [0, 1] but only on 
other parts of the real axis. For instance, from the above discussion it follows that along the 
negative real axis F{z) will have a discontinuity 

AF{x) = F{x + ie) - F{x - ie) = F^^^ (x + ie) - F{x - ie) 

= 2in / w{t)dt, X E negative real axis, (15) 
Jo 



12 



which means that F{z) has indeed a cut along the negative real axis. This definition of the 
'fundamental' Riemann sheet of F{z) coincides with that for the simple example of Fw=i{z) 
discussed in the previous subsection. 

(ii) Alternatively, we could have required that the function F{z) should have no cuts along 
the negative real axis. This amounts to redefining its 'fundamental' Riemann sheet: we start 
again from our zq = Xq — ie with xq < 0, ask that F^^\zq) = F(zo) but then require that the 
values of F^^\z) above the real axis should merge, for Re 2; < 0, with those below the axis. We, 
therefore, define F^'^^z) for Re z < as 

F^^\z) = [\og{t - z)w{t)dt. (16) 
Jo 

This definition which initially has been given only for Re ^ < may then be extended to the 
whole complex plane cut along the real segment [0, 1] (and further, along the positive real axis up 
to infinity; but, as mentioned at the beginning of this subsection, in order to keep the discussion 
simple we shall not consider whet happens beyond the end point z = 1). According to the 
definition (ii), the function F^^\z) will have no cuts along the real negative axis, but in return, 
will have different values along the upper and the lower lip of the real segment [0, 1]. In the lower 
half plane F^^\z) coincides with the function F{z) corresponding to the previous definition (i) 
of the fundamental Riemann sheet, but not any longer for Im z > where F{z) was identical to 
F^^\z). However, since by the construction of the function F^^\z) [see the discussion preceding 
Eq. (ig)] we had 

F^^\x - ie) = F^^\x + ie) + C(e), for e > and < x < 1, (17) 

it follows that F'^^^z) represents now the analytic continuation of F^^\z) on the next Riemann 
sheet — call it sheet (1) — when one crosses the segment [0, 1] in the upward direction. 

On the other hand, starting from the values F*^") [x + ie) from the upper lip of the segment 
[0, 1] and continuing them downwards^ one gets the function F^~^)(z): 

F^~^\x - ie) = F^^\x + ie) + 0{e), for e > and < x < 1 (18) 

which is obtained by deforming the integration contour in the lower half plane. Hence, F^~^\z) 
will have the form 

F^-^\zo) = -2i7i r w{t)dt + \og{t - zo)w{t)dt (19) 
Jo Jo 
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where, again, the logarithm has the same determination as in Eqs. (|T^ and (16). Further, we 
see that the jump of F^'^^ across the cut [0, 1] is 

AF^°\x) = -227r r w{t)dt, for < x < 1. 
Jo 

The functions F^^'^\ F^^\ F^^^ hving on the Riemann sheets (—1), (0), (1) represent in fact an 
unique analytic function the various branches F^~^\z), F^^\z), F^^\z) being nothing 

but its values on a cutting up of the initial Riemann manifold along some arbitrarily given cuts. 

If the function w is holomorphic in some neighbourhood of the origin it admits there an 
expansion of the form 



w{t) = oo + ait + a2t^ + 



and so, from Eq. ([ISD we find 

AF(x) 



2m 



I "1 2 I ^2 3 , 



A function which has precisely the same jump along the negative axis is given by 

z \ogz 



«1 «2 2 I 

aQ + —z + —z + 



(20) 



(21) 



f22) 



which provides a mathematically correct and at the same time extremely simple description of 
the singularity of F{z) near the point z = {]. 

Before closing this section we shall discuss briefly the structure of the singularities at the 
origin of the free terms (72 (-z), g?,{.z), . . . , gj{z) of the iterated equations (p^Ol). In contrast to 
g{z), these functions are written as integrals of the form (|l^) over weights w{t) which are no 
longer holomorphic but contain the singularities of g{z), g2{z), gj_i{z) respectively. By 
straightforward integration we find that the general term of the iterated function gj{z) is 

0™log'=(z), with /c = 0,l,...,j and m > A;. (23) 

At this stage we may wonder whether the left hand cuts of the powers of log(t) appearing 
under the integrals will not interfere with the above analytic continuation process. This does not 
occur since the real part of t becomes negative only along the loop in Fig. 5 around the emerging 
part of the cut of log(t — z) which is always at the same side of the negative real axis, i.e. always 
in the same 'H-ivr zone' of the function log(t). This means that 

A [ \og{t - z) \oi{t)di\ = 2m f t"(log \t\ + mfdt, z G negative real axis. (24) 
Jo J Jo 
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V. SINGULARITIES OF THE SOLUTION OF THE FULL EQUATION 

As mentioned in the Introduction, the specific difficulties of our problem are twofold. We 
have first to handle moving cuts; this question has been largely discussed in the previous section 
devoted to the free term. Second — and this is probably the main difference with respect to the 
classical papers on pinch and end point singularities — we will have to cope with the fact that 
the singularities of the function under the integral are a priori unknown, this function being 
the solution of the integral equation whose analytic properties we are trying to establish. In 
this section we shall address this second problem by solving it first for the eigenfunctions which 
are the natural building blocks of the solution, with the hope that their analytic properties 
(together with those of the free term) will be transmitted to the solution itself. Of course, this 
is not at all obvious since we will deal with infinite series and so new singularities may creep in. 
Therefore before we embark in Section V B on the study of the Riemann sheet structure and 
the asymptotic expansion of the eigenfunctions, we will first make sure in subsection V A that 
the analytic properties of the eigenfunctions do exhaust the holomorphic characteristics of the 
solution. This is probably not entirely pedagogical, but reffects fairly well the way in which our 
work progressed. We shall have often to refer in Section V A to various analytic properties of 
the eigenfunctions which will be proved only later in Section V B. 

This type of analysis presented here is not restricted to this particular integral equation, but 
we hope that it provides a working example of how one could proceed in problems involving other 
kinds of weakly singular kernels. To this end we tried to avoid as much as possible any special 
properties of the logarithmic kernel — for instance the fact that its null-space ker K is empty — and 
show how we can proceed in the general case. At a first reading or if not particularly interested 
in these mathematical proofs but only in the practical aspects of the asymptotic expansion, the 
reader may read only Section V A 1, skip the remainder of the Section V A and pass directly to 
Section V B. 

A. The absence of summation singularities 

We shall proceed in a number of stages: 

In Section V A 1 we discuss the role of eigenfunction expansions in describing the generic 
singularities of the solution. In both Sections V A 1 and V A 2 we recall some facts from the 
theory of integral equations and we identify a class of functions {ip} which can be expanded in 
terms of the eigenfunctions. We discuss also the possible appearance of additional singularities 
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due to problems of convergence of infinite sums of functions. 

The convergence proofs are provided in two separate subsections. In Section V A 3 we 
discuss the continuity of these ■^/'-functions on the real segment [0, 1] and we show that their 
eigenfunction expansions converge uniformly there. However, in order to be able to consider the 
holomorphic properties of the solution f{z) we need a number of results in the complex z-plane. 
These are derived in Section V A4 where, in particular, we prove that no additional singularities 
appear as a result of the summation of the series. The convergence of the asymptotic parts of 
the eigenfunctions is discussed in V A 5. 

1. Generic singularities and eigenfunction expansions 

When discussing the possible singularities of the solution of an integral equation we are 
faced with an apparent paradox. Independently of the specific form of the integral kernel, we may 
always proceed as a numerical analyst usually does when checking the correctness of computer 
code: start in the right hand side of Eq. @ with some 'nice' function f{t) which has no 
singularities, integrate it with the kernel and then choose the free term g{x) to recover the initial 
function /(x). So, irrespective of the (integrable) singularities of the kernel, the solution /(x) 
might be a polynomial, a simple trigonometric function or anything else. One may feel that this 
type of solution is quite exceptional but this example is enough to show that one cannot speak 
about 'the general singularity' of the solution of an integral equation with a given singular kernel. 
However one is fully entitled to ask oneself what may happen in the non-exceptional situations, 
i.e. in the generic case. 

To this end it is enlightening to express the solution of the integral equation for our 
logarithmic kernel in terms of the eigenfunctions Un{x) of K{x,t), defined by 

Unix) = Xn K(x,t)un(t)dt, With. K (x , t) = \og It — x\ , (25) 
Jo 

as the series 

oo ^ 

fix) = gix) + J2 \ 79nUnix), (26) 

n=0 - ^ 

derived below in Section V A 2, where 

9n = f gix)unix)dx. (27) 
Jo 

From expansion ( pB]) it is obvious that in the generic case when small changes in the form of 
the function gix) and hence in {gn} are allowed, both the singularities of the free term gix) 
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and of the eigenfunctions Un{x) will appear in the solution f{x) since they will no longer cancel 
identically. 

2. The functions '^(x) 

In what follows we shall use systematically the notation iIj{x) for the functions from the 
range Ran K of the integral operator. The properties of these functions are used in the derivation 
of expansion (U^ which plays a central role for the analytic properties of / as a superposition 



of those of g and of the u„. However we should exercise great care at this point since additional 
singularities may creep in. We should have in mind the case of the sequence of functions 1, 1 + 
z,l + z + z"^ , . . . . In the open unit disk this sequence tends to the function — z) which 
has a pole at z = 1, whereas all the functions in the sequence are holomorphic in an arbitrary 
large disk. Later in Section V A 4 we will use a theorem of Vitali and/or of Morera to prove 
that no additional singularities appear in a neighbourhood of the origin. The theorem of Vitali, 
for instance, is partially based on the requirement that the elements of the sequence of partial 
sums should be uniformly bounded, which is clearly not valid in the counter-example with the 
geometric series if 1^1 > 1. Hence we have to make sure that in our case all the requirements of 
these theorems are fulfilled. 



In order to obtain Eq. (|26D we first multiply the integral equation (§) by an eigenfunction 
Un{x), integrate over the variable x and use the symmetry of the kernel to get 

fn = 9n + ^/n, (28) 

where the coefficients /„ are defined from f{x) by integrals similar to those in Eq. (P7|). 

From Eq. (^) we find /„ = \ngn/{^n ~ A), but we should avoid expressing the solution 
f{x) as a sum J2 fnUn{x) since the latter might not converge pointwise and/or the eigenfunctions 
{un} might not represent a complete system of orthonormal functions. In the special case of the 
logarithmic kernel (^) it happens (see Appendix) that the do represent such a complete 
orthonormal system, but for an arbitrary kernel K, ker K is not empty and so they do not. 

Many of the textbook theorems concerning expansions of the type (p6|) rely on the continuity 
of the kernel in the square [0, 1] x [0, 1]. Since this is certainly not the case for our logarithmic 
kernel, some supplementary work is necessary to adapt the proofs to our specific conditions. In 
the simple cases where the kernel is continuous one usually takes advantage of this fact to prove 
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that for any square integrable — even singular — function f{t), the function 



ij(x)= [\(x,t)ip(t)dt (29) 
Jo 

from Ran K is (a) continuous for x G [0, 1], and (b) expressible in the form of an uniformly 
convergent expansion 

oo 

^i^) = XI ^""«(^)' (30) 

n=0 

oo 

n=0 



where ipn and are defined by 

i^n= I 'lp{t)Un{t)dt, (Pn= I (p{t)Un{t)dt . 

Jo Jo 

Here the relation which connects the coefficients ipn and cpn can be obtained by multiplying Eq. 
by Un{x), integrating and using the symmetry of the kernel 

= V'n/An. (31) 

The solution f{x) itself does not in general have a representation of the ^/'-kind (|29|), but it is 
clear from the integral equation (P) that the difference f{x) — g{x) is a true ^/^-kind function. 
Hence it can be expanded as the sum J2{fn — fl'n)w„(x), which in turn, using the relation (pSD 



between /„ and gn, yields the representation ( p6D for the solution of the integral equation in 
terms of the free term g{x) and of the eigenfunctions u„(x). 

In Section V A3 we shall prove that the properties (a) and (b) and hence the expandibility 
of f{x) — g{x) are also valid in the case of the logarithmic kernel. Before showing that let us 
notice that if ker K is empty as is the case — see the Appendix — with the logarithmic kernel, or 
if (fi is chosen from the orthogonal complement ker"*" K of ker K, we also have 



= 0. (32) 

L2 



n=0 

However, in contrast with what happens with the -^-kind functions, Eq. (^) represents only 
a convergence in the mean, i.e. <^{x) does not have, in general, an expansion of the form (|30D 
which converges uniformly. For the study of the analytic properties of the solution we shall 
need finer properties than those offered by L^-space arguments, since, for instance, holomorphy 
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and uniform convergence of the partial sums are essential for the Morera theorem to be used in 
Section V A A. 



3. Continuity of xl){x) in the logarithmic case and the uniform convergence of the 
il)^^\x) on the segment [0,1] 



This subsection deals with the properties of the ■j/'-kind functions (|29D on the segment [0, 1] 
for the logarithmic kernel (^). The arguments are quite similar to those which are used in the 
case of the continuous kernels, but we shall discuss them briefly here as a preparation for the 
next subsection and as well as to make this paper self contained. 

Our logarithmic kernel becomes infinite each time the integration variable t equals x. How- 
ever, the continuity (a) of il){x) on the interval [0, 1] (including at its end points) can be proved in 
a straightforward manner using the Schwarz inequality. Indeed, for any function ||v5||l2 < M 
and for any points x and x + 5 belonging to the (closed) segment [0, 1], we have 



(x + 5) - i){x)\^ < f\\og \t-{x + 6)\- log \t - x\]'^dt X 
Jo 



(33) 



where the integral Jq \ip{t)f dt has been replaced by the bound on the L^-norm. If 6 is small 
enough, it can be shown that the integral on the right hand side of (|33|) can be made smaller 
than any given nonzero e VM^, which proves the continuity of all the ip-kind functions of the 
form (p9|). Since the eigenfunctions m„(x) by their very definition (p5[) are also functions of the 
ijj-kind, we have hence implicitly proved their continuity on the segment [0, 1], including at the 
end points a; = and x = 1. 

We shall now show, (b), that for x G [0, 1] the finite sums 



(34) 



n=0 



converge uniformly to the function ip{x) defined in Eq. (pOf) . 

To this end we shall show first that the functions tjj^^\x) converge uniformly to some 
(continuous) function ip^°°\x): a similar proof may be used then for the uniform convergence of 
the analytic extensions ^^^^^z) which will be introduced in the next subsection. From Eq. (pTD 
and from the definition ( pSf ) of the eigenfunctions {u„(x)} we have 



E^IY^-Xnl K{x,t)Un{t)dt 
An 



K{x,t)Y,]^\(PnUn{t)dt 
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which, using the Schwarz inequahty and the fact that the basis {un{t)} is orthonormal, yields: 

EiXl^nUn{x)\^ < £ \K{x,t)\^dt ■ |ESl<^nM„(t)f rft = \K{x,t)\Ut ■ Y.]X\^l- (35) 

For each fixed value of x, < a; < 1, the kernel K{x, t) = log |t — x\ regarded as a function of t 
is in L^[0, 1], and so, 



K^(x,t)dt < < oo. 

JO 



Since the function ip is also in L^, Y^'jXiVn tends to zero for VA; as j increases. Hence the right 
hand side of (^) can be made arbitrarily small irrespective of the value of x. This means 
that the sequence tjj^^\x) converges uniformly to some limit ^^°°\x). Now from the continuity 
of the eigenfunctions which we have proved above, it follows that the finite combina- 

tions {^^^\x)} are continuous. Since the {^^^\x)} converge uniformly, the limit ip^'^\x) is also 
continuous. 

On the other hand, by projecting the kernel onto the basis {«„} we obtain 'the coefficients' 
Un{x)/\n. Bessel's inequality then ensures that any sum over the m^(x)/A^ is bounded: 

E^< / K'ix,t)dt{<oo). (36) 

Integrating with respect to x we see that X)!/-^^ must also converge and hence the |A„| must 
tend to infinity with n. This fact will help us to prove that ^(°°)(x) is in fact identical to the 
function ipi^x) defined in Eq. (^). 

Indeed, since both functions il){x) and il)^^\x) are continuous, it is enough to show that 
the L^-norm of their difference is zero. To this end we notice that the difference between il){x) 
and the functions ijj^^^x) can be written as 

ij{x) - i)^^^ {x) = (x, t)^{t)dt (37) 

Jo 

where K^^'^^\x,t) is the truncated kernel 

K^^^'\x,t) = Kix,t) - ^ (38) 

n=0 ^" 

where we have supposed that the eigenvalues have been labelled according their increasing 
moduli: |Ao| < |Ai| < • ■ ■. Unlike the procedure followed before, we shall not try to use the 
Schwarz inequality to prove directly the pointwise convergence of the ip^^\x), but we shall go 
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instead through L^-space arguments. Since the first j + 1 eigenfunctions Mo(x), Ui{x), Uj{x) 
are all in the null space ker K^-'"^^-' of the truncated kernel (|38|) , it follows that its eigenvalue with 
smallest absolute value is Aj+i. If we denote by K2~^^\x, y) the iterated truncated kernel 

Kt'\x,y) = I K^'^'\x,t)K(^+'\t,y)dt, (39) 



its smallest eigenvalue /i will be A|_,_^. However for any symmetric Hilbert-Schmidt kernel )C{x, y) 
we have 



sup 

\M\=M 



IC{x,y)ip{x)ip{y)dxdy 



where /i is the smallest eigenvalue of JC. Hence, taking the integral over the square of the modulus 
of the left hand side of Eq. (|37D we obtain 

M2 



L2 



^2^^^ {x, y)ip{x)ip{y)dxdy 



< (40) 



which means that \\ip — ip^^^Wi'^ < M/\\j^i\. Since l/|Aj+i| tends to zero as j increases, so does 
11^ _ ^(i) 11^2. Now, from the uniform convergence H^/''--'^ — 'ip^'^m ~^ which has proved above 
[see Eq. (|35D], we have immediately also the convergence 

From the triangle inequality we have 

||^_^(oo)||^, < 11^2 + ||^(^') -^(°°) 11^2, (42) 

where from Eqs. (^Op and (^Tp we see that the left hand side, which is independent of j, can 
be made arbitrarily small by a suitable choice of j in the right hand side. This means that 
11^ — 7/''^°^^ II = which, since both ip^x) and il:'^^\x) are continuous, proves that the two 
functions are identical. 

4. The functions ^^^\z) and the theorems of Morera and Vitali 

In order to be able to study the nature of the singularities of the functions of ■i/'-kind, we 
will have to step off the real axis into the complex z-plane. We shall be particularly interested 
in the complex plane singularities and the asymptotic expansions near 2; = of eigenfunction 
expansions of the form (^), related to the solution f{x) of the integral equation. To this end we 
shall need to know some analytic properties of the eigenfunctions {un}- These will be derived 
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in Section V B where, similarly to the function F^^\z) [see Eq. (0) from Section IV], we shall 
define the analytic functions 

U^^\z) ^=^f A„ log(t - z)un{t)dt, n = 1, 2, 3, . . . . (43) 
Jo 

Although these functions do not yet represent the analytic continuation of the eigenfunctions 
Un{x) which are defined on the real segment [0, 1], they are closely related to them. This specific 
choice is based, as for F^^\z), on the fact that for real negative z we have log(t — z) = log |t — 
since the integration variable t on the right hand side of Eq. (^) is always between and 1. 
In contrast to the function log |t — 2;| (which is not holomorphic because of the modulus), the 
function log(t — z) can be extended analytically in the whole (cut) complex 2;-plane. 

As will become apparent in Section V B, the analytic functions Uj^'^ (z), U^^^ (z), U^^^ (z), 
. . . which are the analogues of the functions F(o)(z), F(i)(z), F^-^\z), . . ., defined in Eqs. (|16D, 
(14) and (19), have an infinite Riemann sheet structure; the superscript in parantheses indicates 
the Riemann sheet under consideration. From the definition of the Riemann sheets and the 
continuity properties across the cut, we have, similar to Eqs.(|l^-(|1^), for any x real between 
and 1: 

Ull'\x + ie) = U^^-'^\x-i£) + 0{e), £ > 0, xG[0, 1]. (44) 

As will be shown in Section V B the eigenfunctions Un{x) defined on the interval [0, 1] can be 
written as simple combinations of the boundary values of f/^'^-*(z), U^}\z) and U^~^\z) on the 
upper and lower lips of the cut: 

Un{x) = ]^{Ul^\x + ie) + U^^\x + ie)), e\0, < a; < 1, 

Un{x) = ^{Ui°\x - le) + Ut^\x - ie)), e\0, < x < 1. (45) 

Hence, in analogy with the function ■ip(x) and the finite sums tp^^\x) defined in the previous 
subsection for x G [0, 1], we shall define now for z in some given region Q of the complex plane, 
the holomorphic functions 

<i/iz) [\ogit-zMt)dt, ^ G L'[0,1], (46) 
Jo 

¥^\z) ± ^„t/l°)(.) ^ ± ^Ui'Hz). (47) 

n=0 n=0 ^" 

Since we are mainly interested in the behaviour of the solution near the origin, it is sufficient to 
take the region Q to be the (open) unit disk cut along the real segment [0, 1] (see Fig. 6), but 
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any other (cut) disk of radius R is also acceptable. The holomorphy of '^{z) in f2 follows from 
the theorem of Morera which states that iff ^^'^^ the function '^{z) is continuous in the open set 
Vt and iff 

/ ^{z)dz = (48) 

JdA 

along the border of every closed triangle A G Q, (i.e. along any reasonable regular closed 
contour), then the function "^{z) is holomorphic in Q. It is clear that these two conditions are 
met by any function having the representation (^61) for z E Q. The continuity of "^{z) can be 
established in a similar way to that of the functions iIj{x) [see Eq. (^3|) ], while the vanishing of the 
integrals ( P5| ) follows — after the interchange of the integral over t G [0, 1] and integrals in the z- 
complex plane — from the holomorphy of log{t—z) as a function of z E Q. The theorem of Fubini^^ 
permits this interchange of the integration order since the function F{z,t) = log(t — z)(f{t) is in 
L^[dA X [0, 1] ] we are interested in. The holomorphy of the functions "^^^^z) is an immediate 
consequence of the fact that they are finite sums of holomorphic functions. 








Im z 
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FIG. 6. The open set f]. 

We denote by \E'°°(2;) the limit of the sequence {\E'*--'^(2;)} wherever this limit exists. We will 
show that there are no new singularities which enter the region as a result of the summation 
process, i.e. that the limit \E'°°(z) is holomorphic in Q. This is an important point: remind the 
counterexample with the geometric series discussed at the beginning of the Section V A 2. 

In studying the analytic properties of the function '^'^{z), the crucial property is again 
the uniform convergence — to be proved below — of the sequence {^(^■)(2)}. This might seem 
surprising since on the real line there exist examples of sequences of infinitely different iable 
functions which converge uniformly to functions which are nowhere different iable. However, in 
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the complex plane the uniform convergence of the sequences can be used in conjunction with 
the theorem of Morera to prove the holomorphy of their limits. Indeed, as a first consequence of 
the uniform convergence of the sequence {"if^^^z)} one obtains the continuity of its limit '^°°{z). 
Secondly, the identity 

/ ^'^{z)dz = (49) 

JdA 

follows from the vanishing of the corresponding integrals over the holomorpic functions ^^^^^z) 
and from the fact that, because of uniform convergence, the integration and the limiting processes 
can be interchanged. 

To prove the uniform convergence of {^^^\z)} for z E we shall proceed similarly to the 
method used on the interval [0, 1]. Since the sums (^7D are finite, they commute with the integral 
(|43|) from the definition of U^^^ (z) and so 

¥^)(z)^j2Y^U!^\^)= C\0g{t-z)j2iPnUn{t)dt- (50) 
n=0 ''^ n=0 

Using arguments similar to those which led to Eq. (^5]), we find 

'< j\^og{t-z)Ut.Y.]tlvl (51) 

where, for all 2; G fi, the integral over the logarithm is bounded while the sum over the coefficients 
tends to zero as j becomes large. This proves the uniform convergence of the sequence 
|v|/0')(2)} in Vl and hence, by the theorem of Morera, that the function '^^{z) has no singularities 
in the region VL. 

These results can also be proved using a theorem of Vitali^^ which states that if: (a) the 
functions ^'^^'^{z) are holomorphic for 2; G fi, (b) the sequence converges uniformly on some subset 
Tj oi VL which has an accumulation point inside Vt and (c) the functions ^!^^\z) are uniformly 
hounded in f2, then the functions ^!^^\z) tend uniformly towards a function ^!^{z) which is 
holomorphic in VL. Note that the subset S may be the segment just above the real interval [0, 1] 
where the uniform convergence has been proved in Section V A?). 

Using Vitali's theorem we can easily find the regions where \l/°°(2;) is holomorphic, by 
looking at the sets where the ^I'^'-^^z) are bounded. In this way we can verify that unlike the 
functions U^\z) which can be continued on higher Riemann sheets, the limit \l/°°(2;) generally 
does not exist there. This is so because the uniform boundness condition (c) is no longer fulfilled 
outside the first Riemann sheet (unless the coefficients ipn decrease exponentially quickly). The 
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reason is the existence of the factor exp(iA„(z — t)) in the higher Riemann sheet continuations 
Ul.^ \z) of the functions U^\z) — see Eq. (165| ) from Section V B 4 — which grows exponentially 
with A„ if {z — t) is complex. 

Finally we shall show that, similar to on the real segment, the limit '^'^(z) coincides 
inside Q with the function "^{z) defined in Eq. (|46|) . This is a direct consequence of the fact 
that the set ker K is empty in the case of the logarithmic kernel. Indeed, replacing Uj^\z) in 
Eq. (|i7D by its definition Eq. (^61) and the Schwarz inequality give 



< 



z 

I \\ogit - z)\^ dt (52) 

L2 





where the right hand side tends to zero when j grows. [For kernels other than the logarithmic 
one with non empty null space, we can define appropriate ^-functions so that the corresponding 
99-functions are contained in ker^K. The simplest way to do this is to begin with a Neumann 
series but stop after few iterations so that the new (^{x) should belong itself to Ran K.] 

5. Sums of asymptotic expressions 

One of the goals of Section V B is to derive asymptotic expressions valid for 2; — > for the 
analytic extensions U^"* {z) of the eigenfunctions: 

U^^\z) = U^:iM + U^:}^^{z). (53) 

The remainder U!;^]^^{z) behaves like 0{\z\'^~^) where k is some given positive integer and 77 > 
but otherwise arbitrary small. However, the coefficients of the asymptotic terms U'fl^y{z) contain 
some (fixed) positive power of the eigenvalue A„, depending on the value of the exponent k. Since 
the {A„} tend to infinity with ra, we should choose carefully an appropriate definition for the 
{^asy} order to secure their convergence. 



The simplest way to solve this problem is to use the iterated integral equations ([T0|) dis- 
cussed in Section III. The eigenfunctions of the iterated kernels Kr{x,t) (|12D are identical with 
those of the initial one, the only change being that the eigenvalues are now (A„)''. This introduces 
a beneficial factor l/X^ in the coefficients ■j/'r-,™ of the new functions 

iprix) = Kr{x,t)^{t)dt, (54) 

which are now 
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FIG. 7. The closure O of $7 drawn in three dimensions in order 
to emphasize that the interval [0, 1] is included twice. 

This ensures the separate convergence of the series 

= E AM%yi^) = t f?f^Ssy(^) (56) 
n=0 n=0 ■^n 

and 

j j ^ 

n=0 n=0 

As a result the asymptotic expansion of the solution of the integral equation will contain terms of 
the form log'^(2;), /c < r, m > /c, as do the iterated free term [see Eq. p3|)] and the asymptotic 
terms Uf^^^{^z). 

The initial range of validity of the asymptotic expressions derived above is the cut open disk 
and so does not yet extend on the real segment [0, 1]. We are of course interested to show the 
correctness of these asymptotic series also on some small real interval < x < xq. As has been 
shown above, the theorem of Vitali fails to work on higher Riemann sheets since the functions 
U''^\z) and lJ^^^[^z) are no longer uniformly bounded and so the sequence {^^^^[^zy\ does not 
converge any more there. However it is interesting and important for what follows to note that the 
sequence {^^^^{z)] related to our integral equation with logarithmic kernel converges uniformly 
also on the closure i7 of 1), i.e. the function \E'°°(z) is well defined and continuous (because of the 
uniform convergence) on the real interval [0, 1], both when approached from above and below. 
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We have used the symbol to emphasize the fact that in all this discussion the open set has 
to be considered as an open subset of the whole Riemann manifold of the solution of the integral 
equation rather than of the 2;-complex plane (whereas normally the closure Q coincides with the 
unit disk \z\ < 1). The set Q contains the interval [0, 1] twice (see Fig. 7) corresponding to the 
fact that the function "^{z) has different limits when z approaches the interval [0, 1] from above 
or below. 

To prove that the sequence {^^^\z)} converges uniformly also on the interval [0,1] it 
is sufficient to note that when z G [0,1], |log(t — z)| is equal to |log|t — z|| if t > z, or to 
I log \ t — z\ ± inl if t < z. Here the sign of in depends on whether z approaches the real axis from 
above or from below. However in both cases the integral from Eq. (^) is bounded and hence, 
the right hand side of (^) can be made as small as one wishes by taking j to be sufficiently 
large. 

It has been already shown in the previous subsection that limit '$°°{z) coincides with the 
holomorphic function "^{z) throughout the open set Q. The limit \E'°^(2;) does not exist beyond 
the two real segments [0, 1] from the border of Q but is continuous up to and on them, because 
of the uniform convergence of the sequence (2)}. The function (z) is hence defined, by 
continuity, in an unambiguous way on the two real intervals on the boundary of Q. It coincides 
there with "^{z), as everywhere else in the closed set Q. 

In this way we have shown that the asymptotic series \E'!^y(z) obtained using the asymptotic 
expansions ofaj{'\z) are valid in a neighbourhood of the origin in Q, and therefore on the (two) 
real intervals [0, xq]. This means that the remainder "^^^{z) of the asymptotic series is bounded 
by 0{\z\''''^) on the real interval [0, a;o] as well as in the open set Q. 

B. Continuation of the eigenfunctions to the complex plane 

In this subsection we shall study the singularities of the eigenfunctions. To this end, similar 
to the analytic extension F^^^ of the free term from Section IV we shall introduce in Section V B 2 
the analytic extensions Uj^'^ of the eigenfunctions Un{x). 

In trying to continue U^^ analytically on higher Riemann sheets, i.e. in trying to construct 
the function U^^ as we did with the free term in Section IV, we face a specific difficulty related 
to the fact that the eigenfunctions m„(x) as they stand, are defined only on the segment [0, 1]. 
This means that we could no longer 'hook' the integration contour, as we did in Section IV where 
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FIG. 8. The eigenfunctions uo{x) and ui{x). 

the weight w(t), being analytic, was well defined not only on the real segment but also on the 
various complex plane deformations of the initial integration path. 

This point will be solved in Section V B 3 where the eigenfunctions Un{x) will be expressed 
as linear combinations of C/^°^ and its Riemann sheet continuations. Another consequence of 
this fact will be the Volterra integral equation which relates U^^ to U^'^ or vice versa (Section 
V B 4). This Volterra integral equation can be solved effectively, providing explicit expressions 
for f/^^^ or U^~^^ in terms of Uj^\ Finally, in Section V B 5 this integral equation is used in order 
to derive the asymptotic series which describe the singularities of the eigenfunctions around the 
origin. 

1. Eigenfunctions of the logarithmic kernel 

It has been proved in the Section V A 3 that, in spite of the logarithmic singularity in the 
integrand, the eigenfunctions {un} defined by 

Un{x) =\n I log \t - X\ Un{t)dt (57) 

■/o 

are continuous functions. The graphs of some of these eigenfunctions are depicted in Fig. 8 and 
Fig. 9. Although these eigenfunctions are finite and look well behaved at the end points and 
1 of the fundamental domain (they are there continuous), their derivatives are there infinite. It 
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FIG. 9. The eigenfunctions ^2(2;), u^ix) and 1*4 (x). 

is the aim of the following subsections to give a full analytic description of the singularities of 
Un{x) at these end points. 

Notice that initially, in order to compute the eigenfunctions Un{x) it has been sufficient to 
take the range of the variable x to be the same as that of the integration variable t. Usually to 
compute the function from the left hand side of an integral equation for values of the variable x 
outside this initial range, we take the so-called Nystr0m continuation, which consists simply of 
substituting the new values of x in the right hand side of the equation. However this is possible 
only if the kernel is holomorphic in some given region and so this method is not applicable in 
our case because of the modulus function inside the logarithm which spoils holomorphy. Thus if 
in Eq. ( [57| ) we take x to be negative, the left hand side is not the analytic continuation 

of the function which appears on the right hand side. 

2. The functions U^^z) and their analytic continuations U!;^'^\z) 

The function which was defined in this naive way for negative values of x, may be 

used to introduce a new function Ul^\z) which is analytic and, although different from Un{x) 
when X G [0, 1], is intimately related to it. Noticing that for a; < and V t G [0, 1] we have 
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log \t — x\ = log(t — x), we shall define U^\z) by 

f/f (^)^=^'A„ /'log(t-^K(t)rft, n= 1,2,3,.... (58) 
Jo 

Taking the cut of log(Z) = log(t — z) along the real negative axis of the complex Z = t — z-plane, 
f/(f)(z) is a complex function of real type, Uj^'\z) = Un\z), having, as we shall show later, a 
branch point at the origin and a cut lying on the positive real axis. Since we are interested in 
the singularity of the function near the origin we shall consider the detailed behaviour of U^^^ (z) 
only in the neighbourhood of ^ = 0. A similar analysis can be performed at the point z = 1. 

The superscript {^^^ in the present case) of U^'^^z) denotes the Riemann sheet under con- 
sideration. Our notation will be such that, by crossing the cut of Ul^'^ between z = and z = 1 
in an anti-clockwise way, the value of the superscript increases by 1: see, in Fig. 10, the full 
line a^^^-P^^^ from the sheet (0), which is continued by the dashed-line a^^^ lying on sheet (1). 
Conversely, crossing the cut between and 1 in the clockwise direction, we pass from the full 
line from the sheet (0) to the dotted-line P^~^^ from the sheet (—1). 
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FIG. 10. Higher Riemann sheet continuations of Un ■ 

We will want to continue the function U^^{z) defined by the right hand side integral from 
Eq. (pSD, starting from some point z = Zq {zq < 0), along the full-line path and the dashed- 
line a^^\ back to zq, obtaining in this way the value W^^zq) of f/„ (z) on the next Riemann 
sheet. But in so doing we will need to deform the integration contour on the right hand side of 
Eq. (|58D into the upper half complex t— plane as we did in the case of the free term in Section 
IV. We are faced, however, with the difficulty that Un{t), as it stands, is defined only on the real 
segment < t < 1. We therefore need to express Un{t) as the value of some analytic function 
on the upper edge of the cut. Here we will be interested only in the upper lip of the cut, since 
moving in an anti-clockwise direction, the analytic continuation path will hook and deform the 
integration contour into the upper half of the complex t-plane (see the curve C from Fig. 5). 
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3. Expressing Un{t), < i < 1, in terms of U^\t + is) and U^\t + ie) 

Our first concern will be to rewrite the functions appearing under the integral sign of Eq. 
([58| ) as combinations of analytic functions. To this end we again use the definition ( [58| ) and start 



from some point z = Zq, zq < 0. We first continue Ul^\z) along a path lying in the lower half 
z-plane to a point Z- = zi — ie, < zi < I, e > 0, below the segment [0, 1] (along the path P^^^ 
from Fig. 10). If the cut of log(Z) = log(t — z) is taken to run along the negative real Z-axis, 
we have 

f log \ t — zi \ + in if t < zi, 
\im\og{t ~ {zi - ie)) = \ 

[log\t-zi\ lft>Zi. 

We now split the integral along the segment [0, 1] into one between and zi (where, of course, 
t < zi and so the integral runs along the upper lip of the cut of the logarithm), and a second 
one, between zi and 1, where the integration points stay apart from the cut of the logarithm. In 
this way we get 

U^^\zi-ie) = Xn / log|t — Zi I Un{t)dt + niXn / Un{t)dt + A„ / log |t — -^1 1 Un(t)dt 
JO Jo Jzi 

= Un{Zi) + iriXn / Un(t)dt (59) 

Jo 

where the last equality follows simply from the fact that the first and the third terms combine 
to give exactly the right hand side of Eq. ( |57| ) . 

Similarly we can make an analytic continuation from z = Zq < (along the path a^'^^ from 
Fig. 10) to the point z+ = zi + ie above the cut, to obtain: 

U^^\zi+te) = Xn / log|t — Zi\Un(t)dt — niXn Un{t)dt + Xn \og \t — Zi\ Un(t)dt 
Jo Jo Jzi 

= Un{Zi) - TXiXn / Un{t)dt. (60) 

Jo 

The integral over M„(t) can be eliminated between Eqs. ( ^9|) and (^) by taking the sum of the 
right hand sides, so that 

un{zr) = l[Ui^Hz,+^e) + Ui^\z,-^e)] + 0{e) 

= l[U^:\z,+^e) + Ul'\z,+^e)] + Oie). (61) 

Here the last line of Eq. (^Tj) follows from the fact that the values of U^^ of the function Un 
on the next Riemann sheet on the upper lip of the cut, merge, by definition (see also Eq. (p!7|)), 
with those of f/^°^ below the cut: 

[/i°)(^i - = U^n\zi + + 0{e) , < ^1 < 1. 
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4. A Volterra equation for U^^^z) 

We shall now be able to deform the integration contour in the complex t— plane in the same 
way as we did in section IV. Recalling that by the analytic continuation of Ul!^\z) along the path 
^(0) _ from Fig. 10 one obtains the function U^\z), we find 



t/y(.„)^A„/^.og(.-,-„) ^-'°'(')+^i''w 



dt (62) 



where the integration contour C is shown in Fig. 5 and where we have replaced n„(t) under the 
integrand by its holomorphic expression (|6lD . The dashed line represents the cut of log{t — zq) in 
the complex t— plane, where, again, the +27r-2one and —in -zone mean that the logarithm differs 
there by +iTT and —in with respect to its mean value across the cut. In the region Re t < 0, 
both [/(^o) and are holomorphic. So — recollect the discussion from section IV which led to 
Eq. ([llD — the integral over the Re t < half-plane part of the contour C yields 

-z7rA„ [\u^^\t) + Ul'\t)]dt, 

J zo 

while the rest of the integral, between t = and t = 1, is identical with that from definition (|58D 
of the function U^\zq). So, Eq. (|62D can be rewritten in the form 

U^r^Hz) = U<^\z) + z7rA„ /'[f/W(i) + U^Ht)]dt (63) 



which is valid for any z in the cut complex plane since it involves only analytic expressions. 
Equation ( |63|) can be regarded as a Volterra integral equation for U^^ (z) if the function f/^°-* (z) 
is known, or equally, as a Volterra integral equation for f/^°^(z) if U!^\z) were known. 

This equation can be solved by differentiation. We find immediately 

or, rearranging the terms, 

dz dz 
This can be rewritten in the form 
d 



dz 



d 



ex.p{+i'KXnz)— exp{—i'rrXnz)U^\z) = exp{—i7rXnz)— exp{+i7TXnz)U^^\ 



dz 



z 
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Now, from Eq. (||) it is obvious that U^^\0) = U^^KO). This is an initial condition which 
permits us to write the solution of the Volterra equation as 

U^n\^) = U^n\^) + 2^^A„ r exp[znXr,{z - t)]U^^\t)dt. (65) 







Similarly we can step backwards and express U!^^ with respect to t/^^-*, or U!^ with respect to 

^ n 

Ui^^\z) = U'^^z) - 2z7rA„ T exp[-27rA„(z - t)]U'^\t)dt. (66) 





5. Asymptotic expansion 

To begin with, we shall suppose that each of the functions U^\z) admits an expansion 
around the origin of the form 

U^:\^) = U^^U^) + U^:iJz), 

= E ^kz' + E E b^,z"^ ■ log^-z) + U^^lJz) (67) 

k=0 m=l k=l 

where U^l^^ is 0{\z\^~^^~^). The choice of these series may look very restrictive, but once we 
have shown their consistency, their uniqueness follows from the uniqueness of the solution of a 



linear integral equation with a Hilbert-Schmidt kernel. In Eq. (|67|) the sum Ylk=o '^^z comes 
from the holomorphic part of U^^ around 2 = 0. Taking as before the cut of log(Z) = log(— 2;) 
along the negative real Z-axis, the right hand side of Eq. (|67|) will be holomorphic throughout 
the domain Q which has a cut running along the positive real axis (see Fig. 6). From the 
definition (EB) we see that U^^^z) is a real-analytic function [i.e. U^\z) = Un\z)] and so the 
coefficients (fc = 0, 1, 2 . . .) and bmk k = 1,2, . . .) have to be real. The coefficients bmk are 
then determined in a recursive manner from the coefficients of the regular part. 

In order to find the coefficients bmk we first notice that the analytic continuation of U^}^y 
across the cut yields 

<Ly(^) = E ^kz' + E E bn.kz"' ■ [log(-^) + 2^^]'- (68) 

k=0 m=l k=l 



Inserting the expressions (|67|) and (|68|) in Eq. (^) we may then compare the various terms 
appearing in the left and right hand sides. This comparison can be done for any z in the open 
set Q, for instance for negative real z, in order to be away from the cut of log(— z). Taking the 
limit z ^ in both sides we remark that 

b,k = 0, k = 2,3,... , (69) 



33 




-0.0016 I ^ ^ ^ ^ ^ ^ ^ ^ ^ 1 

.001 .002 .003 .004 .005 .006 .007 .008 .009 .01 



X 

FIG. 11. One parameter asymptotic fit oIuq^x). 



as a consequence of the fact that the differences of terms of the type z[log{—z) + 2i7T]^ and 
z\og''{—z) from the left hand side of Eq. (6^) yield, by differentiation, terms which tend to 



infinity and which are not compensated by similar terms from the right hand side. Looking at 
the constant terms we get 

bu = A„ao. (70) 
Similarly, by differentiating with respect to z both sides of Eq. (Bl) we obtain 



Jnik 



if A; > m + 1, (71) 



while the coefficients bmk with k < m can be expressed iteratively by means of the coefficients 
Ofc. So, for instance, we find 

621 = , 622 = (72) 

We tried to check numerically the accuracy of this asymptotic expansion in regions close to 
the singularity. To this end we took two terms from the holomorphic part of U^^^y{z), and 
expressed the corresponding bmk coefficients in terms of the first two coefficients Oq and Oi. Since 
Oo is determined by the value at x = of the eigenfunctions which, in turn, is determined by 
the normalization condition, we have in fact only one free parameter, the coefficient ai. Using 



34 




uAx) + 



.001 .002 .003 .004 .005 .006 .007 .008 .009 .01 

X 



FIG. 12. One parameter asymptotic fit of ui{x). 



the above expressions for bn, 621 and 622 in terms of Oq (given) and ai (free), we liave cliosen 
Qi to obtain tlie best fit to tlie first few eigenf unctions, in the region [0, 0.01] near the origin 
where we expected the asymptotic expansion to hold. When we plot the function Un{x) together 
with our asymptotic expansion the two curves appear identical on the interval [0,0.01]. In 
order to show the slight difference in these functions we have defined a new function ■u„(x) = 
Unix) — Un{0) — {un{xi) — M„(0))x/xi which has its end points Un{0) and ^^(O.Ol) at the same 
height and so permits the use of a much enlarged y-axis scale. The corresponding plots for 
the two first eigenfunctions are presented in Figs. 11 and 12. Although as it has been already 
stressed, these curves are just one parameter fits and that moreover wc have restricted ourselves 
only to the first terms in the asymptotic expression, the agreement between these asymptotic 
expressions (the full lines) with the computer calculated points of the eigenfunctions (the crosses) 
is really excellent. 

VI. RESUME AND CONCLUSIONS 

As indicated in the Introduction, the numerical calculations related to the solution of the 
inverse problem for EIT are seriously hampered by the high number of mesh points necessary to 
take into account the sharp peaks of the current density near the edge of the electrodes. Since 
these peaks seemed to be intrinsic objects describable by only a small numbers of parameters we 
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have investigated the details of their analytic structure. To this end we have studied the Riemann 
sheet structure of the eigenfunctions of the dominant singular integral equation relating to the 
solution of the mixed boundary problem for the potential, and derived asymptotic expressions 
both for the eigenfunctions and for the solution of the integral equation. These asymptotic ex- 
pressions provide very simple paramctrisations for the anomalous thresholds, whose effectiveness 
can be judged from the Figs. 11 and 12. 

The paper is constructed as follows: In the Introduction and in Section II the mathematics 
of the EIT modeUing is discussed while the corresponding weakly singular integral equation is 
given in Section III. In Section IV we have discussed the effect on the singularities of the free 
term of the moving cut of the logarithm which 'hooks' the integration contour. The generic 
singularities of the solution of the integral equation are then described as a superposition of 
those of the free term and of the eigenfunctions which were derived in Section V B. Section 
V A contains the proof that although we deal with infinite series, no new singularities appear 
in the neighbourhood of the origin, i.e. the singularities of the solution are really those of the 
eigenfunctions and of the free term. We also determine the limits of the domain where the 
eigenfunction sum converges unconditionally and show how one can extend the validity of the 
asymptotic series also on the real segment [0, 1]. 

We hope that the discussion of the analytic properties of the eigenfunctions and solution of 
this quite special, logarithmic singular equation, will provide a working example which might be 
useful also for the study of the singularities of the solution of other weakly singular integral equa- 
tions. The discontinuities across the cuts will certainly be different, but the general discussion 
will probably be fairly similar. 
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APPENDIX: COMPLETENESS OF THE BASIS 

Following a proof given by G. Aubcrson/^ we shall show in what follows that ker K, the 
null space of the logarithmic kernel, is empty and so the eigenfunctions {Mn}n=o,i,2... do span the 
whole Hilbert space of the functions on [0, 1] . 

Suppose that ker K ^ {0}, i.e. that there exists at least one non zero L^-function such 

that 

/ log \x - y\ 4){y)dy = 0. (Al) 



Defining the function 



we have 



v{x)^^ r^dy, (A2) 
Jo 



dvi Xi 

a) — — — (f){x) almost everywhere, 

b) v(0) = 0. 

We can easily prove that the function v{x) satisfies a Holder condition of index 1/2 

\v{xi) - v{x2)\ < A\xi - X2\^^^ , for Vxi,X2 e (0, 1) , 

where ^4 is a positive constant. 

By integrating the left hand side of Eq. (Al) by parts we find 

/■I /"i v(y) 

/ log\x-y\4>{y)dy = v{l)log{l-x)-V dy 

Jo Jo y — X 

so that from Eq. (Al) we obtain 

r f' I^dy = v{l)log{l-x), (A3) 
Jo y — X 

If we now consider the following function 

F{z) 'M Jz(z - 1) ^^dy for zeD (A4) 
^ Jo y — z 

where D is the complex 2;-plane cut along the segment [0, 1], we can show that: 
{i) F is a holomorphic function in D; 

(in) Im F(x + is) = Jxil - x)V ^ ^^^dy = Jx(l-x)v(l) log(l -x), for x e (0, 1); 

Jo y — X ^ 



(iv) Re F{x + ie) = —7ryx{l — x)v{x), for x e (0, 1). 
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The properties (i) — (Hi) imply that 



^^(1) Vy(l -|/)log (l-Z/) 
y- z 

while from (iv) it follows that for x e (0, 1) we have 



TT Jo y — z Jo 



nJx{l-x)v{x) = ^ ^^—^ -dy - C v{y)dy. (A6) 

^ TT Jo y — X Jo 



v{l)^ Vl/(l -l/)log (l-Z/) 

y - X 

Taking now the limits x\0 and x/l we find 

^(1) 



— log(l - - /' v{y)dy = 0, (A7) 

TT Jo \l y Jo 

^ /' i/^log(l - y)dy - f v(y)dy = 0. (A8) 
TT Jo V 1 — y Jo 



y 

Subtracting (AS) from (A7) we see that v{l) — and hence, from Eq. (A6), 

ttJx{1 — X) •'^ ttJx{1 — X) 

Now, since v{l) is zero, the constant C has to be zero too and so v{x) has to vanish identically 
[it had to be so since neither the right hand side of Eq. (A9) and even less its derivative are L^]. 
This implies that ker K is an empty set and so the eigenfunctions {un} of the logarithmic kernel 
form a complete basis on the segment [0,1]. 
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